First-principles calculation of the piezoelectric tensor d of III-V nitrides 
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We report direct first-principles density-functional calculations of the piezoelectric tensor d 
relating polarization to applied stress for the binary compounds A1N, GaN, and InN. The values of 
d are rather sensitive to the choice of the exchange-correlation functional, and results are presented 
for both the local-density and gradient approximations. A comparison with experiment and with 
values predicted indirectly from the elastic C and piezoconstant *e tensors is also presented. 



77.22.Ej, 77.65.-j 

The piezoelectric tensor d of a polar material relates 
to linear order the induced polarization P to the applied 
stress via 



Pi X/ /, <- rr ' 



(1) 



It is an especially relevant quantity in the field of III-V 
nitride compounds, whose piezoelectric and polarization 
properties are prominent |Q and unusual (2) . As implied 
by its definition, d is relevant to electroacustic applica- 
tions ||, and to the determination of polarization and 
electrostatic fields induced by applied stress or strain in 
devices or epitaxial nanostructures Q]. 

The piezoelectric tensor *e [jl]j^] connecting polarization 
and strain e; via 



13 c j ' 



is related to the d tensor of interest here by 



— dikC, 



(2) 



(3) 



where the Cy are the elastic stiffness constants at con- 
stant electric field. It is thus possible to compute d from 
the knowledge of elastic constants and e-piezoconstants 
(using e.g. the theoretical estimates of Refs. || and M, 
see also below). On the other hand, experiments J7|— [l^] 
directly access the dtj measuring the strain €i caused by 
an applied field E via the converse piezoelectric effect 



(4) 



The dij's in Eqs. |l| and|| are by definition identical Jl3| , 
so that a direct comparison is possible, and especially in- 
teresting and useful. Here we compute directly the dij 
constants in the III-V binary nitrides as derivatives of the 
polarization with respect to stress, calculating the Berry- 
phase polarization of a primitive cell explicitly subject 
to a given external stress. This approach enables us to 
directly compare the calculated values with experiment, 



and with indirect theoretical predictions using Eq. p[ and 
separately calculated C and *e tensors. 

The binary III-V nitrides A1N, GaN, and InN crystal- 
lize in the wurtzite structure, and they possess three inde- 
pendent components of the piezoelectric tensor, namely 
^33j <^3i( = t^32)) an d dis(=d24). To compute these con- 
stants we run a damped Parrinello-Rahman-like dynam- 
ics for each of the compounds considered, with the 
constraint that the system be subject to a given non-zero 
stress. The external stress is applied to the zero-stress 
equilibrium structure, which is identical to that reported 
111 Ref. j§. As a compromise between the contrasting 
needs for sufficiently large stresses to obtain polarizations 
outside the numerical noise, and for sufficiently small 
stresses to avoid large deformations and non-linearity, 
we separately apply stress components a\ = a 2 = ±50 
Kbar in the basal plane, and 03 = ±50 Kbar along the 
singular axis, to determine d%i and (^33; for dig we ap- 
ply a shear stress 05 = ±50 kBar. In the case of InN, 
these stresses and the ensuing deformations may lead to 
metallization because of the very small calculated DFT 
band gap. This must be avoided because macroscopic 
polarization cannot be defined or computed in metallic 
systems. We found that the maximum applied stress had 
to be reduced to 5 Kbar for InN (incidentally, this does 
not causes numerical noise problems in the polarization 
calculations because of the very strong piezoelectric re- 
sponse of InN). After the bulk structure has been op- 
timized for a given stress, the macroscopic polarization 
is computed using the Berry phase technique |Tq ]. The 
stress derivatives of the polarization, i.e. the dy's, are 
computed numerically using a standard two-points for- 
mula. 

As for the technical ingredients, we worked in the 
Density-Functional-theory pseudopotential plane-wave 
framework using the VASP code j[5| for structure op- 
timization and a custom-made code for the Berry phase 
polarization. Ultrasoft potentials |fa] were used for all 
atoms involved , and d semicore states were included in 
the valence for Ga and In. The results reported here 
were obtained with a plane-wave cutoff of 325 eV, and 
an (888) Monkhorst-Pack mesh [|l7j for Brillouin-zone 
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integration; these parameters were found sufficient to 
converge the computed stress for the systems consid- 
ered. As the exchange-correlation functional, we used 
both the Ceperley-Alder |1| Local Density Approxima- 
tion (LDA) in the Perdew-Zunger |L9| parameterization, 
and the Generalized Gradient Approximation (GGA) in 
the variant of Perdew and Wang, known as PW9f pPf . 

In Table | we report all the d^-'s calculated within 
GGA and LDA (the latter are given in square brackets), 
compared with indirect predictions obtained via Eq.[| us- 
ing elastic and e-piezoconstants calculated separately |5| , 
and with the available experimental data [fz| — - The 
calculated results depend considerably on the choice of 
the exchange-correlation functional, the LDA values be- 
ing always larger than those of GGA. One can calculate 
the dij 's using elastic and e-piezoconstant data from e.g. 
Ref [|| and Eq. || to check consistency. These indirect 
prediction are reported for a subset of the dy in Table | 
with the label "indirect", along with the percent devia- 
tions from the directly calculated values. The newly cal- 
culated Cn elastic constants, needed in the evaluation 
of the "indirect" constants, are also listed. The "indi- 
rect" values are always smaller than the directly calcu- 
lated ones, and in most cases they are within ^10% of the 
latter. This suggests that the range of stresses considered 
here is in the linear regime, and proves the (approximate) 
consistency of the data with Eq. |[ 

The available experimental values are as frequently 
above as below or within the GGA-LDA range; the devia- 
tion are in the range of ±4-30 %. From the comparison, it 
cannot be decided whether the LDA or GGA description 
is the most reliable. An exception is the shear constant 
e?i5, which is described appreciably better by the LDA 
for GaN, a similar difference occurring also for InN, for 
which no experiment is available. This difference is pre- 
sumably attributable to an incipient failure of GGA in 
describing accurately the shear elastic response of GaN 
and InN. There are no GGA calculations available for the 
latter; in turn, the LDA values of the shear constants by 
Wright Q is in reasonable agreement with experiments. 

Concerning the level of agreement of predicted and 
measured values, several points should be noted. First, 
all experiments were performed on constrained epitaxial 
samples. The epitaxial constraint affects the elastic and 
piezo response of the epilayer, and as a consequence two 
types of constants are generally reported 0, clamped 
(the one actually measured) and free-standing. The lat- 
ter (with which we compare our calculations), are ob- 
tained from the former via a combination of compliance 
constants (the components of the inverse of C). This 
adds some uncertainty to the comparison, given the large 
spread of the elastic constants values available. The max- 
imum spread for GaN, for example, is as large as 1 pm/V 
using different sets of elastic compliance constants. 

Second, if the substrate on which the piezoelectric sam- 



ple is grown is itself piezoelectric, the total deformation 
measured at the surface of the sample is likely to be 
influenced by the distortion of the substrate. For ex- 
ample, the two GaN samples used in Ref. are re- 
spectively single-crystal (0001 )-oriented GaN grown on 
SiC (0001), a piezoelectric material, and polycristallinc 
(0002)-oriented GaN grown on Si (100); the measured 
clamped constants are 2.8 and 2.0 pm/V in the two cases, 
which translates in free-standing bulk constants of about 
3.7 (see Table and 2.6 pm/V respectively. Assum- 
ing the quality of the samples is comparable, and since 
polycrystallinity in wurtzites is mostly c-axial (vertical 
domains), a polarizable substrate may be argued to lead 
to a ~ 35% apparent increase in piezoelectric response. 
The lesser material quality of the polycrystalline sample 
may be invoked to dispense with it altogether; neverthe- 
less, the argument serves to illustrate the point that the 
typical error bar expected in these experiments may eas- 
ily increase to a few ten percent. Similar considerations 
may apply also to the measurements of Ref. Q where a 
GaN/AlN/Si structure was used. 

Third, a clear result of our calculations is that the ra- 
tio R=— d^/2d^,i is never equal to 2 (it is in fact usually 
larger) . This is to be attributed to the structural and re- 
sponse non-ideality of wurtzite nitrides. In Ref. R, d\% 
was obtained from ^33 assuming that R=2. Our results 
seems to invalidate this procedure, hence the measured 
values of d.13. 

In summary, we have obtained the piezoelectric ten- 
sor d of III-V binary nitrides from first-principles LDA 
and GGA density-functional force-stress and polarization 
calculations. The results are consistent with previously 
calculated e-piezoelectric and elastic tensors. Agreement 
with experiment is only moderately satisfactory (devi- 
ations range typically between 4 and 30%). This may 
be attributed to experimental uncertainties such as the 
translation from clamped to free-standing measured val- 
ues, and the response of polarizable substrates, as well as 
to theoretical uncertainties in the elastic response within 
GGA vs LDA. 
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istry of Research via a Cofin99 project, the INFM Parallel 
Computing Initiative, and a Cagliari University research 
fellowship. 
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TABLE I. Stress-piezoelectric tensor d in III-V binary nitrides, in pm/V. The values listed are: directly calculated within 
GGA and LDA (given in square brackets); calculated indirectly via separately calculated elastic and e-piezo tensors (labeled 
"indirect"; the deviations from the direct calculations are also listed); and experimental results (labeled "exp"). The exper- 
imental cfai constant was obtained Q as ~dss/2 and is not listed here. The experimental value of d-^ was obtained for a 
clamped sample, and the free-standing value reported here was inferred using the elastic response of the material. The Cn 
elastic constant needed in the evaluation of rf lndlrGCt are also listed, in GPa for both GGA and LDA (in square brackets). 



Material 


A1N 


GaN 


InN 


Tmdircct 
"13 


-2.1 [-2.6] 
-2.0, -5% [-2.4, -8%] 


-1.4 [-1.5] 
-1.2, -15% [-1.4, -3%] 


-3.5 [-4.4] 
-3.1, -12% [-3.8, -14%] 



d33 


5.4 [6.4] 


2.7 [2.7] 


7.6 [8.4] 


jindirect 
"33 


5.0, -8% [6.1, -5%] 


2.4, -12% [2.6, -4%] 


6.1, -20% [7.5, -12%] 


,cxp(a) 
"33 


5.6 


3.7 




,cxp(b) 
"33 


5.1 


3.1 




,cxp(c) 
"33 

,cxp(d) 
"33 


6.72 


2.6 




,cxp(o) 
"33 


5.53 






dl5 


2.9 (3.4) 


1.8 (3.3) 


5.5 (8.5) 


,cxp(f) 
"15 


3.6 


3.1 




Cn 


506 [545] 


414 [473] 


266 [314] 



a) See Ref. \ 7\ 

b) See Ref. 

c) See Ref. 

d) See Ref. 

e) See Ref. 

f) See Ref. 



correction for clamping included. 
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